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Abstract. Physical, structural, and electromagnetic properties and interrelating 
processes in sea ice are used to develop a composite model for polarimetric 
backscattering signatures of sea ice. Physical properties of sea ice constituents such as 
ice, brine, air, and salt are presented in terms of their effects on electromagnetic wave 
interactions. Sea ice structure and geometry of scatterers are related to wave 
propagation, attenuation, and scattering. Temperature and salinity, which are 
determining factors for the thermodynamic phase distribution in sea ice, are 
consistently used to derive both effective permittivities and polarimetric scattering 
coefficients. Polarimetric signatures of sea ice depend on crystal sizes and brine 
volumes, which are affected by ice growth rates. Desalination by brine expulsion, 
drainage, or other mechanisms modifies wave penetration and scattering. Sea ice 
signatures are further complicated by surface conditions such as rough interfaces, 
hummocks, snow cover, brine skim, or slush layer. Based on the same set of 
geophysical parameters characterizing sea ice, a composite model is developed to 
calculate effective permittivities and backscattering covariance matrices at microwave 
frequencies for interpretation of sea ice polarimetric signatures. 


1. Introduction 

Sea ice is important to global climate because of its role in 
the mass balance and heat transfer processes between the 
ocean and the atmosphere [Maykut, 1978, 1982; Wettlaufer, 
1991). Global monitoring of sea ice with polarimetric remote 
sensing has drawn considerable interests [Rignot and Drink- 
water , 1994; Cavalieri et al. , 1991]. To interpret sea ice 
signatures in polarimetric remote sensing data, it is neces- 
sary to relate electromagnetic scattering effects of sea ice to 
its physical and structural properties. This relationship pro- 
vides insights into sea ice properties that influence electro- 
magnetic scattering mechanisms and determine polarimetric 
signatures of sea ice. 

Here, the purpose is to relate sea ice characteristics and 
processes to active polarimetric signatures of sea ice through 
the development of a model based on sea ice properties. This 
is presented in a series of two papers: part 1 discusses the 
theoretical model, and part 2 the experimental observations 
[Nghiem et aL, this issue]. In this first paper, sea ice 
properties are discussed and then used to develop a com- 
posite model for understanding polarimetric signatures of 
sea ice including volume and surface scattering mechanisms 
at microwave frequencies. Voluminous polarimetric mea- 
surements of various sea ice types in the Beaufort Sea have 
been archived in the library of data acquired by the Jet 
Propulsion Laboratory multifrequency polarimetric syn- 
thetic aperture radar (SAR). The C band SAR on the first 
European Remote Sensing Satellite has provided backscat- 
tering data of sea ice at small incident angles (20°-26°) [Kwok 
and Cunningham , 1994]. Expeditions and experiments for 
sea ice characterizations spanning several decades have 
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collected numerous physical parameters of sea ice. The 
composite model in this paper provides a link between the 
observed polarimetric signatures to the measured sea ice 
properties. In part II, the model will be used to explain 
experimental observations of sea ice. 

Sea ice has long been a subject of intensive research. 
Extensive results from studies of sea ice and its major 
constituent, ice, have been reported in classical publications 
[Hobbs, 1974; Glen , 1974; Weeks and Ackley , 1982]. The 
complexity of sea ice is attributed to its inhomogeneous 
composition, crystallography, structure, growth process, 
thermodynamic variations, and environmental effects. Inter- 
actions among the physical processes in sea ice result in the 
intricate relationship among physical, structural, and elec- 
tromagnetic properties of sea ice. 

To a limited extent, sea ice properties have been used in 
calculations of conventional backscattering coefficients [Lee 
and Kong , 1985; Tjuatja et al ., 1992] and fully polarimetric 
scattering coefficients [Borgeaud et al 1989; Nghiem et al ., 
1990]. Note that in radar polarimetry, there are other pola- 
rimetric representations [Nghiem et al ., 1990] equivalent to 
the description with the polarimetric scattering coefficients. 
In the above treatments, one or more of the following have 
been assumed: independence of interrelated sea ice param- 
eters, effective isotropy in the sea ice medium, spherical 
scatterers or aligned spheroids, uniform thickness, flat inter- 
faces, no snow or slush or brine layer, and multiyear ice 
without hummocks. These assumptions simplify certain 
properties of sea ice and thus limit the understanding of sea 
ice scattering signatures. 

Affected by inhomogeneities in sea ice, wave propagation, 
attenuation, and polarimetric scattering are interrelated and 
linked to physical and structural properties of sea ice. The 
isotropy assumption contradicts the preferential vertical 
structure in columnar sea ice, which constitutes a significant 
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fraction in the global sea ice mass. Spherical scatterers and 
untilted aligned spheroids cannot account for depolarization 
effects in sea ice polarimetric signatures in first-order scat- 
tering. In the case of tilted aligned spheroids, the calculated 
cross-polarized return has an increasing and then decreasing 
trend as a function of incident angle, while measured data 
show a monotonically decreasing trend. Furthermore, 
shapes of scatterers and their orientations influence wave 
interactions in sea ice even in zeroth-order effects. 

A uniform thickness of an optically thin ice sheet causes 
oscillations in sea ice signatures, while a distribution in 
thickness weakens or even erases these oscillations. Volume 
inhomogeneities in sea ice with flat interfaces give low 
backscattering from sea ice at small incident angles. Rough 
interfaces without volume inhomogeneities underestimate 
scattering magnitudes, depolarization, and decorrelation ef- 
fects at large incident angles. Snow, slush, brine covers, and 
hummocks modify sea ice scattering and need to be consid- 
ered to interpret sea ice polarimetric signatures. 

In the following sections, sea ice physical, structural, and 
electromagnetic properties and their interrelations and vari- 
ations imposed by the physical processes of sea ice growth, 
desalination, and environmental effects are presented. The 
complexity of these sea ice properties is encompassed in the 
composite model. The model includes the thermodynamic 
phase distribution of sea ice constituents, wave speed mod- 
ification, absorption and scattering losses in inhomogeneous 
media, and ellipsoidal inclusions for polarimetric volume 
scattering. Effective anisotropy of sea ice due to the prefer- 
ential vertical alignment of brine inclusions and random 
azimuthal orientations of crystallographic c axes in sea ice is 
considered. Thickness distribution in a sea ice layer, snow 
cover with randomly oriented spheroidal ice grains, brine 
skim and slush layer with high salinity, multiple wave 
interactions with boundaries in anisotropic layered media 
are modeled. Rough surfaces at layer interfaces with wave 
attenuation and differential phase delay through the layered 
media and effects of hummock topography on both volume 
and surface scattering mechanisms are accounted for. Phys- 
ical meanings conveyed in mathematical expressions for sea 
ice signatures are explained. Polarimetric data interpreta- 
tions with the composite model for experimental observa- 
tions of sea ice in the Beaufort sea are presented in part 2 
[Nghiem et a/., this issue]. 

2. Characteristics of Sea Ice 

A background of sea ice characteristics is necessary to 
consider the realistic complexity of sea ice in the modeling 
and to define domains of model input parameters constrained 
by the sea ice physics. This section presents sea ice charac- 
teristics in the context of their interrelationship with electro- 
magnetic properties for the development of a polarimetric 
composite model of sea ice. 

2. 1 . Sea Ice Constituents 

Sea ice is an inhomogeneous medium composed of an ice 
background, brine inclusions, air bubbles, and solid salt. 
Electromagnetic properties of these constituents are charac- 
terized by permeabilities and permittivities, which relate 
material characteristics to electromagnetic fields by consti- 
tutive relations [Kong, 1986|. The sea ice constituents are 
nonmagnetic materials, and their permeability is {4n x 
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10 7 H m^ 1 ). Air has a real permittivity of e 0 (8.8542 x 
10 12 F m~ l ). However, permittivities of other sea ice 
constituents are complex quantities determining wave num- 
bers, with real parts for wave speeds and imaginary parts for 
wave attenuations. 

Ice is a dispersive medium whose complex permittivity 
changes with electromagnetic wave frequency. The imagi- 
nary part of ice permittivity varies by 2 orders of magnitude 
as a function of frequency [Evans, 1965]. For seawater brine 
the absolute value of the permittivity is more than 1 order of 
magnitude larger than that of ice in the microwave frequency 
range. Both real and imaginary parts of brine permittivity 
decrease by several times as frequency increases [Stogryn 
and Desargant, 1985]. Although solid salt volume is usually 
minute compared with the other constituents, the precipita- 
tion process of salt at the eutectic temperature redistributes 
the phase of brine and air inclusions. Mixing of these 
constituents with different phases determines the effective 
permittivity of sea ice. 

With permittivity contrasts to the ice background, inho- 
mogeneities act as scattering sources induced by an incident 
electromagnetic field, giving rise to the source terms in 
inhomogeneous vector wave equations. These current 
sources radiate electromagnetic fields in all directions into 
the upper medium, where the scattered waves are measured 
at all linear polarization combinations. The inhomogeneities 
in sea ice are responsible for effective wave propagation, 
attenuation, and scattering, which are interrelated and linked 
to the properties of the constituents in sea ice. 

2.2. Crystallography and Structure 

In natural geophysical conditions, the ordinary ice poly- 
morph (ice I h ) is predominant [Hobbs, 1974; Glen , 1974; 
Weeks and Ackley , 1982]. Other crystallographic phases of 
ice are formed at very low temperatures or high pressures. 
As illustrated in Figure 1, ice I h has a hexagonal symmetry 
with oxygen atoms packed closely in parallel planes called 
basal planes [Petrenko, 1993]. The crystallographic c axis 
(principal hexagonal axis) is perpendicular to the basal 
planes, and the a axes are parallel to these planes. Natural 
ice is a polycrystalline medium composed of many single 
crystals with different orientations. Beneath a thin transition 
layer, sea ice becomes columnar in structure with c axes 
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Figure 2. Brine pocket sandwiched between ice platelets. 


parallel to within a few degrees of the horizontal plane 
[Weeks and Ackley , 1982]. 

A polycrystal in sea ice consists of small parallel ice plates 
referred as ice platelets. Seawater is trapped in pockets 
between these ice platelets, depicted in Figure 2, during the 
growth of sea ice. As the ice grows and the temperature 
decreases, fresh water in the pockets freezes to the platelet 
walls, leaving liquid in the pockets more saline. These brine 
inclusions have a high permittivity that strongly affects 
electromagnetic properties of sea ice. Horizontal thin sec- 
tions and vertical micrographs of sea ice reveal a substan- 
tially ellipsoidal shape of brine inclusions [Weeks and Ack- 
ley , 1982; Arcane et al., 1986]. Because c axes are parallel to 
the horizontal plane, ice platelets are vertical. Consequently, 
ellipsoidal brine inclusions, sandwiched between platelets, 
are preferentially oriented in the vertical direction. Unless c 
axes are aligned by an underlying sea current [Weeks and 
Gow, 1978], the inclusions are randomly oriented in the 
horizontal plane as shown in Figure 3. 

This orientation distribution of brine inclusions in the 
ordinary ice polymorph renders the effective permittivity of 
sea ice uniaxially anisotropic with the optic axis in the 
vertical direction [Nghiem et al ., 1993b]. The vertical anisot- 
ropy in sea ice effectively supports two characteristic waves, 
ordinary and extraordinary [Kong, 1986], propagating at 
different speeds and attenuation rates. First-year sea ice, 
especially young ice, contains many brine inclusions, and 
their properties are carried in sea ice polarimetric signatures. 
In multiyear sea ice, because of the desalination in the 
growth process (section 2.4), the inhomogeneities are pre- 
dominantly air bubbles, which are more rounded, and thus 
multiyear sea ice becomes more isotropic [Nghiem et al ., 
1993c]. 


2.3. Temperature and Salinity 

Temperature and salinity are two of the most important 
characterization parameters of sea ice. Under stable condi- 
tions, temperature in a sea ice layer is distributed linearly 
from the top interface with cold air to the bottom interface 
with warmer seawater. Many measurements of sea ice 
temperatures show an approximately linear profile in ice 
depth [Cox and Weeks , 1974; Wen et al ., 1989; Gow et al ., 
1990; Eicken , 1992]. In first-year sea ice a C-shaped salinity 
distribution in ice depth is usually observed [Wen et al ., 
1989; Gow et al ., 1990; Eicken , 1992]. Upward brine expul- 
sion and flooding are important mechanisms, causing the 
high top salinity [ Eicken , 1992]. The high bottom salinity is 
obvious because of the interface between ice and saline 
seawater. 

Electromagnetic properties of sea ice are strongly related 
to the temperature and salinity. These two parameters 
together with sea ice density govern the thermodynamic 
phase distribution of sea ice constituents. Fractional vol- 
umes of brine and air inclusions are obtained from equations 
based on phase diagrams [Cox and Weeks , 1983]. Near the 
top of the ice layer the temperature is low and the salinity is 
high (C-shaped salinity profile). Deeper into the middle of 
the ice layer the temperature increases, while the salinity 
decreases. Except near the seawater interface, where the 
brine volume is large, the distributions of temperature and 
salinity effectuate an approximately uniform distribution of 
brine volume in the ice layer [Wen et al ., 1989]. Typically, 
the average brine volume of first-year sea ice in the Beaufort 
Sea is of the order of few percent [Wen et al., 1989; Weeks 
and Ackley, 1982]. 

Furthermore, temperature affects permittivities of individ- 
ual constituents of sea ice. Ice permittivity varies with 
temperature [Tiuri et al., 1984; Mdtzler and Wegmuller , 
1987]. Brine permittivity is a strong function of temperature 
[Stogryn and Desargant , 1985], and seawater permittivity 
depends on both temperature and salinity [Klein and Swift, 
1977]. Temperature also affects sea ice effective permittivity 
by the process of salt expulsion. At the eutectic temperature 
of -21.1°C, a sharp phase change occurs due to the precip- 
itation of sodium chloride dihydrate, NaCI • 2H 2 0 [Weeks 
and Ackley, 1982]. In seawater brine, initial formation of the 
salt appears at -22.9°C [Weeks and Ackley , 1982], and 
consequently, sea ice permittivity drops [Arcone et al., 
1986]. 

2.4. Growth and Desalination 

Ice grows dendritically at large supercoolings due to the 
rapid growth in directions parallel to a axes [Glen, 1974]. 
Columnar ice grows in the form of vertical ice platelets, 
which preferentially align brine inclusions in the vertical 



Figure 3. Horizontal section of sea ice with random azimuthal orientations of c axes. 
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(a) Growth Rate of Sea Ice (March If) 
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(b) Growth of Thin Ice (0-10 cm) 
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Figure 4. Ice growth (a) as a function of thickness 
[Thorndike et al., 1975] and (b) as a function of time. Open 
circles are growth function fitted with data from Thorndike et 
al. [1975] for sea ice, squares are data for ice grown from 
water with 30 %c salinity under a cold plate at -70°C [ Lofgren 
and Weeks, 1969], and triangles are data for ice grown from 
water with 30%c salinity in a cold room at air temperature of 
-22°C during CRRELEX’93. Thickness function h(t) = 
a[ 1 - exp (-ht)] is used for the continuous fitting curves. 


direction, as discussed above. Sea ice growth rate is deter- 
mined by environmental conditions such as air-sea temper- 
ature difference, winds and sea currents, and heat fluxes, 
including latent and radiative heat. 

Thorndike et a!. [1975] have estimated sea ice growth rate 
in polar oceans throughout the year. Figure 4a shows the 
growth rate/(/i) (in centimeters per day) as a function of sea 
ice thickness (in meters) on March 1 1 (winter) with the fitted 
function /(/i) = 12.71 — 0.21 16/r (/? < 50 cm), 2.13 exp 
(-0.017/?) (50 < h < 150 cm), 0.52-8.67 x 10 A h (150 < 
h < 600 cm). From dh/dt = 12.71 - 0.21 16/i (for ice 
thinner than 50 cm), the ice thickness is h(t) = 60[1 - exp 
(-0.21 16/)], where t is the elapsed time from the initial ice 
growth. Thus for thin ice (0-10 cm) the ice growth is 
approximately linear under the first-order expansion of /?(/). 
Figure 4b compares this growth of sea ice with NaCl ice 
grown under a cold plate [Lofgren and Weeks , 1969] and 
saline ice grown in a cold room during the Cold Regions 
Research and Engineering Laboratory Experiment, Septem- 
ber 1993 (CRRELEX’93). For faster growth rates, ice plate- 


let spacing becomes narrower, and the bulk salinity due to 
brine entrapment is higher [ Meese , 1989; Nakawo and 
Sinha , 1981; Lofgren and Weeks , 1969]. This is a determin- 
ing factor in ice crystal sizes and brine volumes. 

As the ice becomes thick, ice growth slows significantly, 
and desalination takes place. There are several desalination 
mechanisms: brine migration caused by thermal gradient, 
expulsion by pressure in brine pockets, drainage by gravity, 
and flushing by surface meltwater [Meese, 1989]. It is found 
that the average salinity in Arctic sea ice is a piecewise linear 
decreasing function of ice thickness [Cox and Weeks, 1974); 
the pronounced change in the slope of the linear salinity 
function of thickness corresponds to the change in the 
dominant desalination mechanism from brine expulsion to 
gravity drainage [Cox and Weeks , 1974]. The desalination 
reduces the lossy brine volume and increases the air volume. 
As a result, the attenuation is decreased, and electromag- 
netic waves can penetrate deeper into the ice. The waves 
encounter more scatterers with weaker local scattering ef- 
fects, since the permittivity contrast of air to ice is less than 
that between brine and ice. This competing effect in the 
scattering from sea ice needs to be accounted for in the 
model to interpret the scattering signature of multiyear sea 
ice, which experiences significant desalination. 

2.5. Surface Effects: Roughness, Snow, Slush, 
and Brine Covers 

Physical properties of sea ice are further complicated by 
surface effects such as roughness, snow, slush, and brine 
covers. Undeformed first-year sea ice is usually smooth with 
some small-scale roughness. Roughness contributes to the 
total backscattering, especially at small incident angles. 
Moreover, multiyear ice undergoes differential melt, which 
results in hummocks with meter-scale lengths. From hum- 
mock profiles of air-snow and snow-ice interfaces of multi- 
year ice in the Beaufort Sea [Cox and Weeks , 1974], large- 
scale correlation functions are obtained and shown in 
Figures 5a and 5b. The hummock topography tilts local sea 
ice surface areas and modulates small-scale roughness. Thus 
hummocks modify both the volume and surface scattering 
properties of sea ice. 

Snow cover on sea ice introduces an additional scattering 
layer and also causes physical changes in sea ice. While the 
sea ice layer is anisotropic, the snow medium is effectively 
isotropic because of random orientations of ice grains in 
snow [Sihvola and Kong , 1988; Nghiem et aL, 1993b]. The 
isotropic snow cover has a masking effect, which lessens the 
anisotropic characteristics in the polarimetric signature of 
columnar sea ice. Furthermore, snow capillary force gives 
rise to brine wicking, especially for thin ice [Drinkwater and 
Crocker , 1988]. Consequently, a slush layer with high salin- 
ity appears between the snow and the ice layer. This slush 
layer has a high permittivity that strongly influences the 
scattering from sea ice. A highly saline surface skim, or a 
thin brine layer, has been observed on the surface of young 
lead ice [Perovich and Richter-Menge, 1994; Drinkwater and 
Crocker , 1988]. The brine layer significantly impacts polari- 
metric signatures of sea ice. 

3. Polarimetric Composite Models 

Based on the properties of sea ice discussed above, a 
composite model for polarimetric signatures of sea ice is 
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(a) Hummock Snow-Ice Interface 



(b) Hummock Air— Snow Interface 



Horizontal Displacement (m) 

Figure 5. Normalized correlation functions of hummock 
profiles for (a) snow-ice interface and (b) air-snow interface. 
Data are derived from hummock profiles of snow-covered 
multiyear ice in the Beaufort Sea [Cox and Weeks , 1974], 
continuous curves are for Gaussian correlation functions, 
and dash-dotted curves are for exponential correlation func- 
tions. 


developed and presented here. Both effective permittivities 
and polarimetric scattering coefficients are derived from the 
same set of sea ice properties. Their interrelationship is 
therefore accounted for and linked to sea ice physical and 
structural properties and processes. 

Sea ice covered by a layer of snow, slush, or brine is 
modeled with a layered configuration, as illustrated in Figure 
6. The upper half is air with permittivity e 0 and permeability 
Ijlq. Region 1 is described as an isotropic inhomogeneous 
medium containing randomly oriented spheroidal scatterers 
with a spatially dependent permittivity £\(r). Sea ice in 
region 2 is an inhomogeneous medium, with permittivity 
e 2 (r)> composed of an ice background with ellipsoidal brine 
inclusions whose orientation distribution renders the me- 
dium effectively anisotropic. For multiyear ice, air bubbles 
are included in the ice layer. Region 3 is the homogeneous 
seawater with permittivity e 3 . Amplitude vectors described 
in Figure 7 for upgoing and downgoing waves will be used to 
derive multiple wave interactions with boundaries in the 
layered media. 

Effective permittivities of inhomogeneous media are ob- 
tained with the strong permittivity fluctuation theory. Cova- 
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Figure 6. Layered configuration of sea ice. The top layer is 
brine, slush, or snow cover; the middle is sea ice; and the 
bottom is the underlying seawater. 


riance matrices, whose elements are polarimetric scattering 
coefficients characterizing polarimetric signatures of sea ice, 
are then calculated with the analytic wave theory, which 
preserves the phase information. Small-scale roughnesses 
and large-scale hummock topography are also considered. 
Wave scattering mechanisms in the inhomogeneous aniso- 
tropic layered media with composite interfaces are illus- 
trated in Figure 8 and modeled in the following sections to 
derive polarimetric signatures of sea ice. 

3.1. Effective Permittivities of Inhomogeneous Media 

The mixing of the constituents in the inhomogeneous sea 
ice determines an effective permittivity, which governs wave 
propagation and attenuation in the medium. In this model the 
dispersion of sea ice not only depends on the dispersive 
permittivity of each constituent but also on the scattering 
effect of the inhomogeneities in addition to the quasi-static 
part. Effective permittivities of sea ice are derived from the 
strong permittivity fluctuation theory [Tsang and Kong , 
1981], which is extended to account for randomly oriented 
spheroidal scatterers [Yueh et al., 1990; Nghiem et al. y 
1993a], and for ellipsoidal scatterers [ Nghiem et al . , 1993b]. 



Figure 7. Amplitude vectors in the layered configuration. 
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Incident Wave a 



Results for effective permittivities are summarized here, and 
their physical meanings are pointed out in terms of sea ice 
properties. 

Consider medium n, with n = 1 for the cover layer such 
as snow and n = 2 for the sea ice layer (see Figure 6). 
Generally, to account for the medium anisotropy, the effec- 
tive permittivity is described by the tensor 

£ neff & ng £<neff* * £ neff 0) 

where / is the unit dyad, e ng is the auxiliary permittivity, 
is the dyadic coefficient, and £„ eff is the effective dyadic 
scatterer [Tsang and Kong , 1981]. 

Mathematically, the first term, e ng , in (1) together with the 
dyadic coefficient S n are determined by the condition of 
secular-term elimination. The secular term originates from 
the singularity corresponding to the coincidence of an ob- 
servation point and a source point in the inhomogeneous 
medium. The results for e ng can be reduced to a generalized 
form of the Polder-van Santen dielectric mixing formula 
[Tsang et al ., 1985]. 

Physically, e ng is the effective permittivity under the 
quasi-static condition. It characterizes wave propagation in 
an effective medium without the scattering effects. From the 
functional forms of e ng [Nghiem et al ., 1993a, b], it is 
observed that 


£ ng ^ ng^f nr* £ nc & ng' (2) 

where subscript c denotes a constituent in the inhomoge- 
neous media. As indicated in (2), e ng is a function of 
fractional volume f nc of the medium constituents, permittiv- 
ity contrast (e m . - and dyadic coefficient element 

ISnl 

In sea ice (n = 2) the constituents are ice, brine, or air as 
discussed in section 2.1. The fractional volumes of these 
constituents follow the thermodynamic phase distribution 
determined by the temperature, salinity, and density of sea 
ice (section 2.3). Constituent permittivities are complex and 
vary with frequency, temperature, and salinity (sections 2.1 


and 2.3). In first-year sea ice, brine inclusions and the ice 
background have a large permittivity contrast (section 2.1), 
which strongly affects wave propagation and attenuation in 
sea ice. In desalinated multiyear sea ice (section 2.4) the 
permittivity contrast between air bubbles and ice is smaller, 
wave speed is larger, absorption is smaller, and wave 
penetration is deeper compared with the case of first-year 
sea ice. 

Dyadic coefficient S n does not depend on the size of the 
inclusions but does depend on ratios of ellipsoidal axial 
lengths and orientation distributions of the nonspherical 
scatterers [Nghiem et al ., 1993b]. Physically, it is the shape 
and orientation that determine the polarizability of the 
inclusions. Thus the ellipsoidal shape and orientation of 
brine inclusions (section 2.2) are manifested in electromag- 
netic properties of sea ice even in the lowest order. 

The second term in (1) accounts for the wave scattering 
effect in sea ice. This term depends on frequency; therefore 
the inhomogeneous sea ice is inherently dispersive (even if 
one assumes frequency independence for the constituent 
permittivities). In the existence of inhomogeneities, wave 
speed is modified and wave attenuation is increased due to 
the scattering loss. Effective dyadic scatterer in the 
second term of (l) is calculated with 


[f/ieffL 



x,y,z 

d<t>nfPn('l'nf> <t> nf) 2 F 
kj 



dP[G ag (k')] k & n€ (k') + [S n ] kl 




(3) 


in which k 0 is the wave number in air, G ng is Green’s 
function, r ( n °^ kim is the variance of scatterer tensor £ related 
to permittivity fluctuations in the inhomogeneous medium 
[Tsang and Kong , 1981; Nghiem et al . , 1993b], <t> n ^ is the 
Fourier transform of the normalized local correlation func- 
tion R n ^ describing the size and shape of the scatterers, and 
Pni&nf' <f>n /) * s the probability density function of orientation 
with Eulerian angles i(t n f and <t> n f relating local scatterer 
coordinates to the global coordinates (x, y, z). Subscript £is 
used to indicate that T^ /m corresponds to The vari- 
ance r<% im has the same functional dependence as in (2) 
and thus the same physical relationship to sea ice properties. 

For orientation distributions, p„( «/*„/, 4> n f ) = sin (^ rt /)/ 
(47 t) describes random orientations of scatterers such as ice 
grains in snow, and p n (ip n f , <£„/) = £(<//„/)/( 2 7 r) depicts brine 
inclusions oriented in the vertical direction (delta 8 function 
of (//„) in columnar ice with random c axes in the horizontal 
plane. The ensemble average over random orientations in 
snow leads to an isotropic effective permittivity. For the 
sea ice the effective permittivity is vertically uniaxial, in 
conformance with the preferential alignment of brine 
inclusions. 

The integration over the product of Green’s function and 
the correlation function depends on inclusion sizes; there- 
fore the size effect appears in the higher order compared 
with the shape and orientation effects. A normalized local 
correlation function is used to model the ellipsoidal brine 
inclusions in sea ice. In the local coordinates (jc\ >*', z’) 
associated with the principal axes of an ellipsoid, the corre- 
lation function is 
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with correlation length l nx >, l ny > , and ( n = 2) related to 

the effective size and shape of the scatterers [Nghiem et al ., 
1993b]. For spheroidal scatterers such as ice grains in the 
snow layer (n = 1 ), the correlation function is reduced from 
(4) with = /„ v . = l„ p •. 

After substituting correlation function (4) into (3), the 
integrations are carried out over the complex integrands for 
to obtain effective permittivity from (1). Expres- 
sions for effective permittivities with spheroidal and ellipsoi- 
dal scatterers are reported by Yueh et al. [1990] and Nghiem 
et al. [1993a, b]. Note that for sea ice ( n - 2) composed of 
ice and brine, fractional volumes of the constituents are 
related by / 2ice = 1 - fi brine- In this case » Nghiem et al. 
[1993b] use the notation with subscript s for the scattering 
brine (f 2s = /Wine and e 2s = £ 2brine) and subscript b for the 
background ice ( f 2 b - 1 " fix and e 2h = ^ 2 ice)- In these 
calculations, appropriate branch cuts and Riemann sheets 
are considered for the involving functions with complex 
variables. The result for ellipsoidal scatterers can be reduced 
to the case of spheroids by equating two of the three 
correlation lengths. For randomly oriented scatterers, the 
effective permittivity becomes isotropic, and is the unit 
dyad / multiplied by a scalar permittivity e neff . 

3.2. Volume Inhomogeneities in Layered Media 

Inhomogeneities in sea ice, besides influencing wave prop- 
agation as presented in section 3.1, give rise to wave 
scattering from sea ice into the upper media. Scatterers such 
as ice grains in the snow layer also contribute to the 
scattered fields. Thus physical and structural properties of 
sea ice are related to both wave propagation and scattering 
with interactions to the cover layer. 

These processes are included in the polarimetric model 
with the analytic wave theory applied to the layered config- 
uration in Figure 6. Scattering sources, due to the permittiv- 
ity contrasts of the constituents in the inhomogeneous me- 
dia, are introduced in wave equations for the layered media. 
Scattered fields are solved from these equations under the 
distorted Born approximation [Tsang and Kong , 1981], with 
effective permittivities and dyadic Green’s functions ac- 
counting for wave interactions in the media. The set of sea 
ice characterization parameters used in the calculation of 
effective permittivities is utilized in the derivation of the 
scattered fields. Wave propagation, attenuation, and scatter- 
ing are therefore consistently linked to sea ice characteris- 
tics. 

Scattered fields from sea ice are measured in the backscat- 
tering direction with a polarimetric scatterometer or radar 
for all combinations in the linear polarization basis. Ensem- 
ble averages of scattered field correlations are used to obtain 
the complete set of polarimetric backscattering coefficients. 
These coefficients constitute a covariance matrix or a Muel- 
ler matrix, characterizing polarimetric scattering properties 
of the layered media [Nghiem et al., 1990]. To derive the 
scattering coefficients, operator L„j klm is defined for conve- 
nience as 


x.y.z 



i.j'k.i.m 



dd„ q„(d„) 


i: 


d'l'nf 



nf Pnitff nfy nf 


n C >\y Gi ’ 



dr n 


fi nfy fi nf ) 


(5) 


where q n is the probability density function of thickness d n , 
p n is the probability density function of orientation, and 
location r n or is in volume V n occupied by region n. The 
orientation distributions are the same as discussed in section 
3.1 for effective permittivity. The thickness distribution 
characterized by q n is governed by a first-order differential 
equation depending on sea ice growth rate, motion, and a 
redistribution function related to ice deformation [ Thorndike 
et al . , 1975]. 

The correlation function C ,,^/m of the scatterer with 
variance T njkim is determined by the following expression: 

C n£jklm{r 1 1* Gj > *1* nfy fi nf 


j dp r nJUn&nfiP) e X P [-'73 • {ft n ~ Ol ,6) 


where <2> n { is the Fourier transform of (4) and £ n j k is the local 
scatterer [Nghiem et al 1993a, b]. 

The operator (5) is applied on elements (G) of mean 
dyadic Green’s functions and mean incident fields ( F ) to 
obtain the correlations of scattered fields under the distorted 
Born approximation as follows 

2 

<£ 0i (r) * E* 0s (r)) = ^ L nikl m l(C {]mj (r, r„))(F nk (r „ )>] 

n = 1 


• [<G 0n „(r- OXFnJKW <7) 

This result includes the physical and structural properties of 
the anisotropic sea ice, wave interaction with the medium 
interfaces, and effects of the cover layer. The results given 
by (7) can be reduced for application to a simpler sea ice 
configuration without the cover layer by setting thickness d { 
to zero. 

The mean dyadic Green’s function (G 0 „> in (7) is for the 
observation point in region 0 and the source point in region 
n. The Green’s function has been derived for an anisotropic 
medium in a two-layer configuration [Lee and Kong , 1983] 
and extended for application to a polarimetric model of 
layered geophysical media [Nghiem et al . , 1990]. The dyadic 
Green’s function in the radiation field can be written as 

e ik " r 

(C 0 „(f, r,)) = - — exp (~ik ■ p s )g„(k p , z s ) (8) 
47T/' 


>i = 1,2 

where e lki ' r /(4nr) is the spherical radiation, k p is the lateral 
wave vector, and p s — xx s + yy s . Dyadic coefficient g, ,(A p , 
z v ), representing multiple wave interactions with medium 
boundaries in layered media [Nghiem et al.. 1990], is defined 


as 
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9„(k p , Zs) = 2 2 £(*0 Z ) 

fl V 

•W nil v(-£ p )Hk™) exp ( ~ik™z s ) 

+ V H ^(-E,)Hk«) exp (-ikgz,)] (9) 

where polarization unit vector jx can be h or — z>, subscript p 
stands for horizontal h or vertical v polarization, and super- 
scripts u and d denote upgoing and downgoing directions, 
respectively. In region 1, unit vector v — (i , -v and 
subscript v - h, v without a minus sign in front of v. Since 
region 1 is isotropic, there is no distinction in the propaga- 
tion of horizontal or vertical waves; thus /:(*" = k\ z = k“ z 
and k ^ = k\f = kf z . In region 2 the medium can be 
anisotropic, and characteristic waves are ordinary o or 
extraordinary e. In this case, unit vector v - o, -e and 
script v - <?, e. There are four types of waves in region 2: 
ordinary downgoing, ordinary upgoing, extraordinary down- 
going, and extraordinary upgoing. Thus the birefringent 
effect (double refraction), wave speed difference, and differ- 
ential attenuation in the anisotropic sea ice medium are 
accounted for in this model. All polarization and wave 
vectors are defined by Nghiem et aL [1990] and U and D will 
be shown below. 

For incident field F 0j - = [ft(k 0zi )E hi + exp 

(/& 0 ,* • r) y the mean fields can be written as 

(F n (r)) = exp (ik pi « p)P n (k ph z) n= 1,2 (10) 


B„ = D„ • Bn 


(13b) 


Elements in U n are U tlfXl/ and in D n are D npin where p and 
yean be h or v except that v — o t e in region 2. Matrix £/ 0 
is defined as the reflection matrix R 0 for region 0, and Z> 3 is 
the transmission matrix t 3 for region 3. From the system of 
matrix equations in (12), downgoing and upgoing amplitude 
vectors are solved in terms of Z? 0 , and the results are then 
compared with (13) to obtain coefficient matrices in the 
dyadic Green’s functions and the polarization vectors. Phys- 
ically, these vectors and matrices describe wave interactions 
with the interfaces including multiple reflections and trans- 
missions in the sea ice and the cover layer (such as snow, 
slush, or brine), the reflection to the upper air, and the 
transmission to the underlying seawater. 

By applying operator L from (5) on the dyadic Green’s 
functions of the layered media and the mean fields as in (7), 
correlations of the scattered fields are obtained. Then pola- 
rimetric backscattering coefficients <J plVK are calculated with 


JU. T VK 


lim 

r—*9c 

A—** 


47 rr 2 

^ E Ti E*. 


(14) 


where subscripts fx , v , r, and k can be h or t\ subscripts i and 
s stand for incident and scattered fields, respectively, r is the 
distance from the radar, and A is the illuminated area. 
Coefficients from Greens’s functions and mean fields are 
combined to arrive at the following expression of the scat- 
tering coefficients in the scattered polarization basis: 


where subscript / indicates the incident wave, p = xx + yy 
is the lateral spatial vector, and C pI = xk xi + yk yi — A 0 ( Jc sin 
0 O( cos <f) 0i + y sin 0 O( - sin <£ 0( ) is the lateral component of 
incident wave vector k 0h 0 Oi is the incident angle, and <f) 0i is 
the azimuthal angle. Polarization vector P n (k pi , z) is ex- 
pressed as follows 

= E pi [Dn P Ak pi )Hk%) exp (ik^z) 

M v 

+ U npv (K pi )HkZ) ^P dk v n “z)] (11) 

where \x = h y v and v — h y v in region 1 or v = o y e in region 
2. As indicated in (11), the incident wave is effectively 
decomposed into the ordinary and extraordinary waves in 
the anisotropic columnar sea ice. 

Coefficients U and D are derived from boundary condi- 
tions. To illustrate the wave interaction processes described 
by U and D, consider amplitude vector A n of upgoing waves 
and B n of downgoing waves in region n = 0, 1, 2, 3. 
Amplitude vectors of waves propagating away and toward 
each interface, as shown in Figure 7, are related with matrix 
equations 


' ' 


B n{n+ 1 ) 

T{ n +\)n 


B n 

B{n + 1) 


Tn(n+ l ) 

B(n + \)n 


^(/i + n 


In region 0, i? 0 is the amplitude vector of the incident wave. 
In region 3 there is no upgoing wave and R 32 and T 32 are not 
needed. Amplitude vectors in different layers can be related 
to the incident vector B 0 by 

A n = U n B 0 (13a) 



x,y\z 

■P n (<knJ- <!>„)) 2 2 

,d 


(15) 

in which «, b — — 1, I for n = 1 and a y b y c\ d — on , od y 
eu , ed for n — 2. 

In (15) the contribution from the isotropic region 1 in = 1 ) 
consists of 16 terms for downgoing and upgoing incident and 
scattered waves. In region 2, four wave types in Green’s 
function and four in the incident fields provide 16 terms, 
which are correlated with 16 terms in the scattered waves, 
giving rise to 256 terms from the sea ice layer, accounting for 
all types of wave interactions in the anisotropic layered 
media. Details of and $ n are given by Nghiem et al. 
[1990, 1993b]. The integrations over Eulerian angles in (15) 
can be carried out numerically. When the low-frequency 
condition is valid, analytical expressions can be derived. For 
a given thickness distribution, the integration over the thick- 
ness probability density function gives the polarimetric 
scattering coefficients. For a uniform thickness, q n becomes 
a Dirac delta function of thickness, and the results are 
readily obtained. 

Constituents of sea ice are nonmagnetic (section 2.1), and 
thus the medium is reciprocal. The cover layer and the 
underlying seawater are also reciprocal. The reciprocity 
reduces the backscattering covariance matrix to a 3 x 3 
complex Hermitian matrix. In the covariance matrix, con- 
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ventional backscattering coefficients <r hh = & hhhh , cr vv = 
<r vvvv , and <r hv = cr hvhv are the diagonal elements. The 
cross-correlation coefficients are a hhhv and cr hvvv . The co- 
polarized ratio is defined as y = o- vv la hh , the cross-polarized 
ratio is e = o- hv /a hh , and the complex correlation coefficient 
between the horizontal and vertical returns is p = 

&hhvv( &hh 

In sea ice where c axes have random orientations in 
azimuthal directions, the reflection symmetry completely 
decorrelates the cross terms ( o- hhhv = a hwv = 0), and the 
covariance matrix for sea ice contains only five independent 
parameters at oblique incident angles [Nghiem et al. , 1993c]. 
In (9) and (1 1), and k™ are different from k%z and *2z in 
both real and imaginary parts in the uniaxial sea ice. Conse- 
quently, the scattering centers of horizontal and vertical 
waves are separated because of different wave speeds and 
attenuations. Then, p calculated from the scattering coeffi- 
cients with (15) has a reduction in magnitude and an increase 
in absolute phase compared with the isotropic case. At the 
normal incident angle, the covariance matrix has two inde- 
pendent parameters, since y = 1, Im p = 0, and e = (1 - 
p)/2 due to the azimuthal symmetry [Nghiem et al., 1992]. 
These relations in polarimetric coefficients of sea ice are 
followed by (15). When c axes are aligned by a sea current, 
sea ice becomes a biaxial medium, which is not treated in 
this paper; the polarization signature of a biaxial sea ice is 
generally characterized by a covariance matrix with nine 
independent parameters. 

Effects of a cover layer with different compositions and 
characteristics such as snow, slush, or brine (section 2.5) on 
polarimetric signatures of sea ice can be seen from the 
results in this section. For snow-covered sea ice, (15) is the 
sum for the scattering from snow (n - 1) and sea ice (rt = 
2). Due to the isotropy, horizontal and vertical returns from 
snow are better correlated. The overall effect increases the 
magnitude of p and decreases the absolute value of the phase 
compared with the case of bare ice. These polarimetric 
characteristics approach the behavior of an inhomogeneous 
isotropic medium, and thus the anisotropy of the sea ice 
underneath is partially masked. 

When sea ice is covered by a high permittivity layer such 
as slush or brine, wave attenuation is increased because of 
the highly lossy saline layer. Furthermore, magnitudes of the 
elements in /? 0 i increase, and the difference between the 
horizontal and vertical copolarized elements an d 

Rq\ vv are larger, especially at larger incident angles. In this 
case, backscattering coefficients decrease and copolarized 
ratio y increases at larger incident angles. Further complex- 
ity is introduced by rough interfaces in sea ice, which are 
considered in the next section. 

3,3, Rough Interfaces in Layered Media 

Natural interfaces in sea ice are rough with various length 
scales (section 2.5). In this section, surface scattering from 
rough interfaces in layered media is presented for small-scale 
roughnesses with small standard deviation heights and 
slopes; large-scale roughnesses such as hummocks are ex- 
amined thereafter. To estimate effects of rough interfaces on 
polarimetric signatures of sea ice, contributions from the 
rough surface scattering are considered with wave interac- 
tions, differential propagation delay, and wave attenuation in 
the anisotropic layered media. 

As depicted in Figure 8, the incident wave is scattered by 


the rough interface between air and snow. The wave contin- 
ues propagating into the snow layer (region 1), where both 
the coherent and incoherent fields are scattered again by ice 
grains and the rough snow-ice interface. Part of the wave is 
transmitted into the ice layer (region 2) and is further 
scattered by brine inclusions in this layer and by the rough 
ice-water boundary. The wave also reflects back toward the 
top interface and is scattered again by the inhomogeneities 
and the rough surfaces. This process of multiple scattering is 
complicated and involves surface scattering from above and 
below the interfaces, multiple-volume scattering in the inho- 
mogeneous layers, higher order interactions between vol- 
ume and surface scattering mechanisms, and wave propaga- 
tion and attenuation in the multilayered anisotropic 
inhomogeneous media. The estimation of rough surface 
effects here is simplified by considering the first-order scat- 
tering of the downgoing incident field and the higher order 
scattering is ignored. 

For the top interface the first-order polarimetric backscat- 
tering covariance matrix seen in region 0 is 


C m — 


■ (°> 

°Q\hhhh 


n- {0) 

a 0\vvhh 


rr (0) 

a 0\hvhv 


n- (0) ~ 

^Qihhvv 


rr (0) 

*■' 0 1 ww 


(16) 


This is added to the total scattering [Lee and Kong , 1985], 
implying the volume and surface scattered fields are not 
correlated. Rough surface backscattering coefficients 
aofhhhh and can b e computed with the Kirchhoff 

approximation, geometric optics approximation, or small 
perturbation method [Tsang et al., 1985]. Under all of these 
first-order approximations, the cross-polarized return 
o-tfll.hv is zero and the correlation term tr^ vv = = 

l/2 - 

For the rough surface at the middle interface between 
snow and ice, wave propagation and attenuation in the 
layered media are considered. The matrix method with wave 
vectors (section 3.2) is employed in the derivation of surface 
scattering with wave interactions in layered media. Instead 
of using (Icpi) and ( k p ) to denote coefficients corresponding 
to incident and scattered waves, subscript i is used for 
incident and s for scattered as shorthand notations. As 
illustrated in Figure 7, incident wave amplitude B x in region 
1 is transmitted from region 0 by matrix D u operating on f? 0 
in region 1; explicitly, B { = D Xi * B Q . Then B x is scattered 
by the rough surface characterized by S 12 , which renders 
scattered wave amplitude A x = S 12 * B x . Scattered wave 
A , goes upward from region 1 into region 0 by transmission 
matrix T ias through the air-snow interface; thus A 0 = 
T XOs • A,. The phase factor and attenuation in region 1 of 
thickness d x (see Figure 6) can be expressed in the back- 
scattering direction as exp (ilc u * r u ) exp (i£ Xs • r Xs ) = exp 
[i2(k pi p + k l { ziR d i)] exp (-2 k? ziI d x ), where p is the lateral 
distance and scripts p, z , R , / , and u stand for the lateral 
component, vertical component, real part, imaginary part, 
and upgoing direction, respectively. Under consideration of 
wave propagation, attenuation, and interactions, surface 
backscattering due to the middle interface observed in region 
0 is described by 


SiT = t,o.v • ^12 • Du exp [i2(k pj p + k'{, JR d ,)] 


•exp (-2* jV,) 


(17) 
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In the anisotropic layered configuration with vertical optic 
axis under the aforementioned approximations in the deri- 
vation of rough surface scattering, horizontal and vertical 
waves are decoupled, and the cross-polarized returns van- 
ish. For the observation in region 0, nonzero elements of the 
covariance matrix C 12 obtained from (17) for the backscat- 
tering from the rough interface between regions 1 and 2 are 

a \2hhhh = (^\0hhs^\hhi^\2hh) 

* (T \0hhsD IhhiS \2hh)* ex P (~^k" zii d l) 


In the vertical uniaxial layered configuration, horizontal 
and vertical components are decoupled. In this case, the 
necessary elements U ohi and U evs are obtained from (19) 
and expressed as 


Uohs ~ 


* - e ak ' :,d 'R [lhh R lohh 

e ilk, “- k ' 1 “ )d 'T 2 i ev T l o vv 
I - e ak '» d 'R nvt Jt l0vv 


(20a) 

(20b) 


— \T \0hhsD lhhi\* ex P ( 4ki zi /d \)<T \ 2 hhhh (18a) 

O’Uv vvv = ^ tOvvs^ I2ui) 

■(T ]Qvvs D ]vvi S nvv )* exp {—4k“ z jjdi) 

= I T 10 vvJ> l vvi\ 2 exp ( -4 k “ ziI d ,) a , 2vvvv ( 1 8b) 

a nhhvv = (T\0hhsD\hhiS]2hh) 

{T\Qvvs^ lvv$ 1 2 vi') exp \ z n d \) 


where k 2zs = k 2zs = ~k 2 t$- Furthermore, wave propagation 
in region 2, whose interface with region 3 locates at the depth 
of d 2 , imposes on the backscattered wave a phase factor and 
attenuation of 

exp (ilc 2oi * r 2i ) exp (i£ 2os * r 2s ) 

= exp 0*2* pr -p) exp [i(k< 2zsR + k$ ziR )(d 2 ~ d x )} 

•exp [ — (k 2zs j + k 2ziI )(d 2 — d])] (21a) 

exp (ilc 2ei • r 2i ) exp (i£ 2es ■ r 2s ) 


- ( T \ohh S £) ihhi)(T \Q VVS D ivvi )* 

* exp (-4k" ziX di)<ri2khw (18c) 

where (7\2hhhH> &i2vvvvi an d &\2hhvv — ( cr \2hhhh ( *]2vvvv) 

are polarimetric backscattering coefficients for the rough 
surface. In the calculation of cr 12/An . K , incident angle 0 U in 
region 1 is determined by 6 U ~ tan -1 (k pi ik lziR ) due to the 
phase matching conditions. Under the Kirchhoff or geomet- 
ric optics approximation, the Fresnel reflection coefficient at 
the interface is evaluated at normal incidence without dis- 
tinguishing horizontal or vertical wave polarization. In the 
small perturbation method the two wave types are treated 
differently. Since region 2 is vertically uniaxial, the reflection 
coefficient is computed with an ordinary wave for horizontal 
polarization and an extraordinary wave for vertical polariza- 
tion. 

Scattering coefficients of the rough surface at the bottom 
interface between ice and water are derived in a similar 
manner. Incident wave amplitude B 0 in region 0 transmits 
into region 2 as fi 2 = D 2i ■ B 0 . This incident wave then 
undergoes the scattering characterized by 5 23 from the 
rough interface between regions 2 and 3. Scattered wave 
A 2 = Sj 3 • B 2 is observed in region 0 as A 0 = U s • A 2 . 
Matrix U s is derived by considering wave transmission and 
reflection in the layered configuration. Wave amplitude 
vector A | is transmitted into region 0 as A 0 = Tj 0s • A,, 
where A j in region 1 is composed of the reflection R l2s * B j 
from the interface between regions 1 and 2 and the transmis- 
sion T 2Xs • A 2 from region 2. Thus A { = R ]2s • + 

T 2U * A 2 , which relates A j to A 2 as A x = (/ - 
B i2j * R\o s )~l * ? 2 \s ^ 2 2 where / is the identity matrix. 

By operating T m on A ] , A 0 is obtained in terms of A 2 , and 
U s is identified as 


U K = T 


10 * 


(I ~ Ri2s *^10 s) 


T 21,* 


U ohs U ovs 
U ehs U evs 


(19) 


where o denotes ordinary wave and e is for extraordinary 
wave. 


= exp (i'2 k piP ) exp [i(k e 2lsR - k 2ziR )(d 2 - d x )] 

• exp [ — (k 2 “ s2 — k 2 d jt )(d 2 — ^1)] (21b) 

for ordinary and extraordinary waves, respectively. 

After wave propagation and multiple interactions with 
boundaries in the anisotropic layered media are considered, 
the covariance matrix observed in region 0 has the following 
nonzero elements in terms of backscattering coefficients 
derived from (20) and (21) for the rough interface between 
regions 2 and 3 

0-23 uu = \UohsDzhoi\ 2 exp [-2 (k° 2ls , + k 0 2zil ){d 2 - </,)] 


•exp [-2(k lzsl + k Xzil )d\]cr 22hhhh (22a) 

^ 2 3 Lw = evs^2vei\ 2 ex P [~2(k 2zs/ ~ k 2ziI )(d 2 — d x )] 

• exp [ — 2 (k XzsI + ki ziI )d]]a 22vvvv (22b) 

a 23hhvv = (Uohs&2hoi)( U evs® 2ved* 

• exp [i{k 2zsR + k 2ziR - k 2zsR + k 2ziR )(d 2 — d |)] 

• exp [-{k 2zsI 4- k 2ziI + k{ zs f - k 2z if)(d 2 ~ d\)] 

•exp [-2 (* Ui/ + k j z ii)d j ] cr 

23hhvv (22c) 

Under the Kirchhoff, geometric optics, or small perturba- 
tion approximations, <r 2)M „ = (a- nhhhh a 23m , m ,)' 12 . In the 
calculation of the scattering from the interface between 
regions 2 and 3, the ordinary wave reflection coefficient is 
used for horizontal polarization, and the extraordinary wave 
for vertical polarization. Incident angles in region 2 for 
ordinary and extraordinary waves are determined by 0 2oi = 
tan -1 (k pi /k 2z ift) and 6 2ei = tan" 1 (~k pi /k 2 d iR ), respec- 
tively. Note that there is an extra phase term in <r$ hhvv (22c) 
due to the difference in propagation speeds of ordinary and 
extraordinary waves. 

Surface scattering dominates the total backscattering at 
small incident angles where the information carried by the 
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scattering from the ice layer underneath is masked out of 
polarimetric signatures. Surface scattering depends on the 
roughness and the permittivity contrast between the upper 
and lower media at the medium interface. For the same 
roughness, a higher permittivity contrast gives higher back- 
scattering. In snow-covered sea ice, the permittivity contrast 
at the air-snow interface, especially for dry snow, is smaller 
than that at the snow-ice interface. In this case, the contri- 
bution from the snow-ice rough surface is more important, 
and the total sea ice signature carries some information 
pertaining to this surface. 

Roughness at a sea ice interface evolves under many 
processes related to environmental conditions such as atmo- 
spheric boundary layer, winds, sea currents, wave actions, 
solar radiation, and heat transport. These processes create 
many types of sea ice with complicated surface conditions, 
compositions, and structures such as composite pancake ice, 
rafted ice, pressure ridges, puddles, and thaw holes, which 
are subjects of future model developments. In the following 
sections, melt hummocks, another complication in polari- 
metric sea ice signatures, are modeled and presented in 
terms of their effects on both volume and surface scattering 
mechanisms in sea ice. 

3.4. Volume Scattering With Hummock Surface 

Multiyear sea ice survives at least one summer season 
during which differential melt creates hummocks on the ice 
surface with a large-scale roughness of the order of a few 
meters [Cox and Weeks , 1974; Weeks and Ackley , 1982]. 
During winter, multiyear ice is covered by a new layer of 
snow over the hummocks. Both the top interface between air 
and snow and the middle between snow and ice have this 
undulating topography. 

Evidently, a local incident angle of the electromagnetic 
wave differs from the global incident angle set by the antenna 
looking direction. Effects of this meter-scale topography on 
the volume scattering from the layered media are accounted 
for by averaging scattering coefficients over the distribution 
of local incident angles. This averaging process implies the 
characteristic length of the hummock topography is large 
enough that the local scattering can be treated indepen- 
dently. The method has been applied to composite models 
for rough surface scattering from small-scale roughness with 
large-scale modulation [ Wright , 1968]. When surface slopes 
of the upper and lower hummock interfaces are highly 
correlated with a positive correlation coefficient, the total 
volume scattering from the layered media can simply be 
averaged over the distribution of the local incident angles at 
the top interface. Under this condition, the two interfaces 
are closely in parallel, and their large-scale characteristic 
lengths are similar. 

Scattering coefficient due to the volume scattering 
under the effects of the hummock topography is deter- 
mined by (t^Xk = / 4>i) dflfj dCl, which is the 

average of the local scattering coefficient , evaluated 
at local incident angle 8 ( and azimuthal angle <f> h over the 
observation solid angle ft. For the unit vector h = (— x\ ~ 
yip + z)i(x 2 + + D ,/2 normal to the local surface and 

backscattering direction k = (-Jr sin 0 O/ cos 4> Qi - y sin 0 ()/ 
sin <f> 0i + z cos 0 O/ ) defined by the global incident angle 0 Of 
and azimuthal angle </> 0 , , the average for c t^Xk can be written 


(T 


(V) 

\L VT K 



dx dip o-'^vrJOi, 4>,)p(x . <P) 


• [I + tan 0 Oi (x cos <t> oi + <p sin <j> w )] 


(23) 


in which p {\ , <p) is the probability density function of slope 
and local incident angle 0/ is 


cos 0 / = k • h — 


cos 0o/ + sin 0 O /(a cos </>o / + <P s * n 

~( 2“7 2 , n 1/2 
KX + <P + 1) 

(24) 


Since the scattering configuration has reflection symmetry 
with respect to all vertical planes, the azimuthal angle (f> oi 
can be set to zero without loss of generality. In this case, the 
limit of integration for ip is from — » to +°° and for x is from 
-cot 0 Oi to +oc. The lower limit -cot 0 O/ is introduced to 
account for the shadowing effect of the large-scale surface 
[Chan and Fung , 1977]. For a Gaussian distribution of slope 
S = 2 l,2 a R il R of the large-scale topography with standard 
deviation height a R and correlation length l R , the slope 
distribution is p(x , <p) = exp [-(a 2 + <p 2 )/(2S 2 )], and 
polarimetric backscattering coefficients under the hummock 
effects are determined by 


r (V) = 


/: 

•/: 


e ~X*H2 S l ) 

dx (1 + X tan 0 O <) — — 


d(p (r'„ TK (0i)e 


-* 2 /(2S~) 


(25) 


When the local surface is tilted, the incident horizontal 
polarization is composed of both transverse electric and 
transverse magnetic components. This can be interpreted as 
a roll of the incident polarization basis by an angle a between 
the incident horizontal vector ( t = z x k t l\z x k t \ and the 
local horizontal polarization ( t = n x kjl\n x k t \ (see Figure 
8) for global vertical unit vector z and incident unit vector 
kj = —k. All scattering coefficients are expressible in 
terms of sin a and cos a [Nghiem et al 1992]. For <£ 0f - = 0, 


sin a = 


(hi x fi t ) * kj = 


2? 

[<P 2 + (sin 0 O / — X cos 0o/) 2 ] l/2 


(26a) 


:os a = f i • hi = 


sin 0 O / - x cos 0( 


f (D 2 + (sin 0n; - y COS 0n;) 2 l 


2i 1/2 


(26b) 


Locally, all cross correlations in <r MrrK are zero due to the 
local reflection symmetry. In the global coordinates, all 
cross correlations in cr^lx are a Is° expected to be nullified 
due to the global reflection symmetry [Nghiem et al ., 1992]. 
This can be demonstrated by considering scattering coeffi- 
cient oj'fc/nn for instance. After the polarization roll and the 
reciprocity are used, 

v’hhhv = sin a cos a[-a hhhh cos 2 a + a vvvv sin 2 a 
+ <r hhmi cos 2 a - a vvhh sin 2 a 
+ 2cr/ n ,/ 1I ,(cos 2 a - sin 2 a)] 


as 


(27) 
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Since sin a is an odd function of <p, while cos a is even, the 
integration over <p renders ajXlv zero when p(\, <p) is an 
even function of <p, which is the case for a nondirectional 
slope distribution of the large-scale topography such as the 
Gaussian distribution. All other cross correlations in d^J TK 
are shown to be zero in the same manner. For the calculation 
of the nonzero cr^J. K , the necessary local polarimetric scat- 
tering coefficients are 

v'hhhh = °hhhh cos 4 a + c r vvvv sin 4 a 


statistically independent and stationary Gaussian processes. 
The Gaussian distributions of probability densities for these 
azimuthally symmetric roughness profiles are prescribed by 

I ( fr\ 

Pr(fr) = - 1/2 exp - —2 (29a) 

(2tt) o, \ 2<T r J 

Pr(/r) = - .1/2 — exp (29b) 

(2tt) o r \ 2 o R 


+ (Re & hhw + 2cr^,, At ,)2 cos 2 a sin 2 a (28a) 

® vvvv ~ ^ vvvv cos 4 a + o hhhh sin 4 « 

+ (Re o hhvv + 2o hvhl ,)2 cos 2 a sin 2 a (28b) 

<r'hhw = o hhm , cos 4 a + o vvhh sin 4 a 


where cr r and cr R are the standard deviations of the small- 
scale and the hummock profiles, respectively. Under these 
conditions, the probability density p{f) of the total rough- 
ness profile is a convolution of p r (f r ) and p R (f R ) [Papoulis, 
1984]: 

p(f)=p(fr+fR) 


+ (<rhhhh + <r„ V w) cos 2 a sin 2 a 
— 4 <r hvhv cos 2 a sin 2 a 


(28c) 



df R p Xf ~ f r) P R^f r) 


v'hvhv = °hvhv( cos 2 a - sin 2 a) 2 

+ (<Thhhh + Vvmm - 2 Re o hhm ) cos 2 a sin 2 a (28d) 

As seen in (28), has a contribution from cr M ^ K with a 
factor of cos 4 a. Usually, a is small and this term is large 
compared with the terms involving sin a. Consider the 
cross-polarized return <r’ hvhv (28d) for multiyear sea ice: the 
term (cr hhhh + cr vvvv - 2Re cr hhvv ) has the same order of 
magnitude as cr hvhv [Nghiem et al ., 1993c]; therefore cr' hvhv 
is dominated by the contribution from the first term in (28d) 
for small angle a in the trigonometric factors. For large-scale 
hummocks with a long correlation length compared with the 
standard deviation height, mainly retains the charac- 
teristics of with some polarization mixture from other 
coefficients. In general, hummock effects on volume scatter- 
ing are to reduce both polarization and incident angle depen- 
dences in sea ice signatures as indicated by (28) and (23). 

3.5. Surface Scattering With Hummock Surface 

Effects of hummocks on backscattering from rough inter- 
faces can be estimated by averaging local scattering coeffi- 
cients over the slope distribution of the large-scale topogra- 
phy. If the small-scale roughness is valid for the small 
perturbation method, scattering from small surface pertur- 
bations is first calculated and then averaged over the large- 
scale slope distribution. Then the results are added to the 
scattering from the large-scale roughness obtained under the 
geometric optics approximation. This approach uses differ- 
ent approximations together and considers small-scale 
roughness being tilted locally over patches of large-scale 
roughness according to the slope distribution of the large 
surface. In contrast, the method used in this paper considers 
the small-scale and large-scale roughnesses at the same time, 
applying KirchhofFs approximation to the total roughness, 
and calculating the total scattering from the composite 
surface altogether. Here, KirchhotFs method, which has 
been used to obtain results for single-scale roughness [Beck- 
mann and Spizzichino , 1963; Tsang et al ., 1985], is extended 
to derive backscatter from the composite rough surface. 

Let the total profile / = f r + f R be the superposition of a 
small-scale profile f r and a hummock profile f R , which are 


(2 7r) V2 (a 2 + a 2 R ) 


exp 


f 


2((r} + cr R ) 


(30) 


which is also a Gaussian distribution whose height variance 
c t 2 = cr } + a R is the sum of the individual variances. 

The joint probability density p(/j ,/ 2 ) for two points at p] 
and p 2 on the surface is [Davenport and Root , 1958] 


P(/i> fi) 


I / A -2C/,/ 2 +/| \ 

2tj-o- 2 (1 - C 2 ) l/2 CXP \ 2o- 2 (l-C 2 ) ) 


(31) 


where C is the normalized correlation function between the 
two points. For the azimuthally symmetric stationary sur- 
face, C depends only on the distance p' = |pj — p 2 |. From the 
total profile / = f r + f R} C is derived as follows: 

o 2 C(p') = (fipXfiPi)) = ojCXp') + o R C R {p') (32) 


To arrive at (32), the independence between f r and f R and 
the following correlation functions have been used 

(fXP\)fXPi)) = o 2 r CXp') = oj exp (-p' 2 // 2 ) (33a) 

</*(P.)/tf(P2)> = o 2 R C R (p') = tr 2 exp (-p' 2 //£) (33b) 


where l r is the correlation length of the small-scale rough- 
ness and l R is that of the hummock topography. As seen in 
(32), the normalized correlation function C is a linear com- 
bination of small-scale correlation function C r and hummock 
correlation function C R . 

For a composite surface, consisting of large-scale hum- 
mocks and small-scale roughness characterized by the above 
descriptions, polarimetric backscattering coefficients a { ^J TK 
are derived with the stationary phase method under Kirch- 
hofF s approximation. The results for nonzero backscattering 
coefficients for the composite rough interface at oblique 
incident angle 0, are 


* 2 |/?ol 

7 T COS 2 8 



NGHIEM ET AL.: POLARIMETRIC SIGNATURES OF SEA ICE MODEL 


13,677 


in which ^ and v can be h or v polarization, k is the wave 
number in the upper medium, |/? 0 | 2 is the Fresnel reflectivity 
between the upper and the lower media at normal incident 
angle, and the quantity D is defined by 

D= f d(f> j dp' p'{exp [ — k dz a 2 {\ — C(p'))] 

Jo Jo 

- exp (~k dz a 2 )} exp (ik dp p ■ p') (35) 

where the subscript d denotes the difference between the 
incident and scattered wave vectors; thus £ d = k ( - lc s = 
pk dp + zk dz in the cylindrical coordinates. In the backseat- 
tering direction, k dp = 2 k 2 sin 0 ( and k dz = 2 k 2 cos 0,. 

Carrying out the integration over azimuthal angle <f> results 
in the Bessel function J^(k dp p) and then substituting the 
correlation function C(p’) from (32) yield 


» j 

,1m ( m + 1 \ 


II 

q 

M 

-e \ 2j a r R a mj A m , 


m= 1 | 

V=' J 



(36) 


where a mj is the binomial coefficient and A mJ is an integral 
over p' determined by 

A mj =2 P° d P ' pVo(^ p p')C r (p') m " a ' ,) C J »(p') U " ,) 

Jo 


roughness is blended together, and a joint probability density 
function (31) of the resultant roughness is used. Here the 
surface scattering is not generated independently first at two 
different scales and then added together. Instead, result (39) 
takes care of the fact that there is roughness at different 
scales in the total profile. Consequently, (39) contains terms 
depending on correlation lengths and standard deviations of 
both small- and large-scale roughness components. The 
above results for polarimetric backscattering coefficients 
from rough surfaces under effects of large-scale hummocks 
analytically reduce to the small-scale surface scattering 
when the large-scale surface profile is removed by letting a R 
- 0. Furthermore, the results approach the geometric optics 
solution in the high-frequency limit. At large incident angles 
the Kirchhoff results for rough surface scattering become 
less accurate; however, the surface contribution already 
becomes small and the volume scattering takes over the total 
sea ice signatures. 

In the layered configuration, wave propagation and multi- 
ple interactions with layer boundaries are taken into account 
as shown in section 3.3, where o- n i n + \)n.TVK is replaced by the 
corresponding cr^*, with a R and l R replaced by cr ri(ri + n/? 
and /„(„+!)/? to denote the parameters between medium n 
and ( n + 1), in the calculations of observed in 

region 0. Since the composite interfaces contain different 
scale lengths, large-scale hummocks and small-scale rough- 
nesses have different responses at different wave frequen- 
cies. 


1 

exp 



(37) 


To obtain the result in (37), (33) has been used for the 
correlation functions and v m j has been defined in terms of 
the correlation lengths as 


m-(J-i) j- 1 

V mj _ ,2 + i2 

l r l R 


(38) 


With the above solution, the result for backscattering coef- 
ficient cr { ppvl > from (34) can be written compactly as 


r (S) _ 


kW 

cos 2 Of 


« m + 1 

E X (k dz a mj ) lm a m jA K 

m = 1 7=1 


(39) 


where a m j = and cr m j are determined by the 

following equations, 


j 

1)]! 


(40) 


a mj ~ a r a R 


(41) 


respectively. 

In the above approach, the hummocky surface of sea ice, 
which has a roughness with a small scale riding on the large 
scale of the hummocks is described as a composite rough 
surface by superimposing the profiles of the small- and 
large-scale components. The resultant probability density 
function of the total profile is a convolution of those of the 
individual roughness components as given by (30). Thus the 


4. Summary and Discussion 

This paper presents a composite model for polarimetric 
signatures of sea ice. The purpose is to relate sea ice 
physical, structural, and electromagnetic properties and 
processes to sea ice polarimetric signatures. The paper 
presents new developments in volume and surface scattering 
in the modeling of sea ice such as thickness distribution, 
effects of polarization mixing, and incident angle variations 
due to hummocks on the scattering from volume inhomoge- 
neities, differential phase delay and attenuation for surface 
scattering in layered anisotropic media, and the mixing of 
small roughness with large hummock topography. Based on 
sea ice characteristics, effective permittivities are derived 
with the strong permittivity fluctuation theory, and polari- 
metric backscattering coefficients are obtained with the 
analytical wave theory under the distorted Born approxima- 
tion. Surface scattering from rough layer interfaces is also 
considered with effects of wave attenuation and phase delay 
in an anisotropic layered configuration for sea ice. 

Sea ice is an inhomogeneous medium composed of various 
dispersive constituents. These constituents are nonmagnetic 
materials preserving the reciprocity in sea ice signatures. In 
the form of the ordinary ice polymorph, natural polycrystal- 
line ice is composed of multiple ice platelets. These platelets 
entrap brine in ellipsoidal pockets with preferential vertical 
alignment. Random orientations of crystallographic c axes in 
azimuthal directions render sea ice uniaxially anisotropic. 
The anisotropy causes an effective birefringence, differential 
phase delay, and attenuation in sea ice. Temperature, salin- 
ity, and density determine the phase distribution of sea ice 
constituents, including salt expulsion and brine solidification 
below the eutectic temperature. From the phase distribution, 
fractional volumes of sea ice constituents are obtained to 
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derive both effective permittivities and polarimetric scatter- 
ing coefficients of sea ice. Thus effective permittivities and 
scattering coefficients are consistently linked to sea ice 
characteristics. 

Dendritic growth of sea ice forms columnar ice where 
crystal size and brine volume depend on the growth rate, 
thereby influencing wave propagation, attenuation, and scat- 
tering. Desalination processes in sea ice increase wave 
penetration while weakening local scattering effects due to a 
low permittivity contrast between air bubbles and ice back- 
ground. Isotropic snow cover can partially mask the aniso- 
tropic characteristics in sea ice signatures. Furthermore, 
snow affects salinity distributions by brine wicking and, 
consequently, changes both effective permittivities and po- 
larimetric coefficients. Slush or brine layer with a large 
permittivity inflicts large reflection and attenuation, de- 
creases scattering contribution from the ice layer, and in- 
creases the copolarized ratio at large incident angles. 

Effects of roughness are important at small incident an- 
gles. Rough surface scattering depends on both roughness 
and permittivity contrast at the interface. The contribution 
from a rough interface beneath the ice layer suffers from 
differential attenuation and phase delay. Hummocks have an 
averaging effect on polarimetric volume scattering, which 
mixes local polarization responses and reduces the depen- 
dence of sea ice signatures on incident angles. Moreover, 
hummocks modulate the small-scale roughness, and the 
response from the composite interface is frequency depen- 
dent. From symmetry properties of sea ice scattering con- 
figurations, the number of independent parameters con- 
tained in covariance matrices for polarimetric signatures of 
sea ice are also determined. All surface scattering models in 
this paper are for first order, which is sufficient for field 
measured parameters used in part 2. In this case, the surface 
scattering itself has negligible depolarization and decorrela- 
tion effects for backscattering. 

Based on physical, structural, and electromagnetic prop- 
erties and their interrelating processes, the composite model 
can be used to interpret polarimetric signatures of various 
ice types. For bare sea ice the results are obtained by setting 
the thickness of the cover layer to zero. For sea ice in a 
frozen lead with a brine skim or slush cover, polarimetric 
signatures are calculated with a layered configuration where 
the ice layer has a thickness distribution and the cover layer 
has a large effective permittivity determined from its highly 
saline composition. For snow-covered first-year sea ice, 
signatures obtained from the model with snow and sea ice 
properties include scattering from snow with spheroidal ice 
grains, sea ice with ellipsoidal brine inclusions, anisotropy 
due to sea ice structure, wave interactions with boundaries, 
and rough interface effects in the anisotropic layered sea ice 
media. For multiyear sea ice, the model accounts for hum- 
mocks on the snow-covered ice with surface tilting effects 
and composite rough interfaces. 

In part 2 of this series, the model will be used to interpret 
measured polarimetric signatures of sea ice in the Beaufort 
Sea. Sea ice characterization data obtained from sea ice 
expeditions and experiments will be utilized. The trends in 
experimental data will be explained, and the utility of 
polarimetric signatures for global sea ice monitoring will be 
discussed. Further than a simple suggestion of using polari- 
metric data as a parametric classification of ice types, the 
scope of these two papers provides certain understanding of 


polarimetric signatures based on sea ice physics. For higher 
order effects in sea ice signatures and more complicated sea 
ice types, properties, and processes, further developments 
in sea ice modeling and measurements are necessary. 
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